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Abstract. It was shown recently by Cordoba, Faraco and Gancedo 
in [T] that the 2D porous media equation admits weak solutions with 
compact support in time. The proof, based on the convex integration 
framework developed for the incompressible Euler equations in [4], uses 
ideas from the theory of laminates, in particular T4 configurations. In 
this note we calculate the explicit relaxation of IPM, thus avoiding Ti 
configurations. We then use this to construct weak solutions to the un- 
stable interface problem (the Muskat problem) , as a byproduct shedding 
new light on the gradient fiow approach introduced by Otto in [14) . 



1. Introduction 

We consider the incompressible porous media equation (IPM) in a 2- 
dimensional bounded domain C M^. The flow is described in Euler ian 
coordinates by a velocity field v{x,t) and a pressure p{x,t) obeying the 
conservation of mass and the conservation of momentum in the form of 
Darcy's law: 

dtp + div (pv) = 0, (IPMl) 
divz; = 0, (IPM2) 
t, + Vp = -(0,p). (IPM3) 

Here we chose xi as the horizontal and X2 as the vertical direction, with 
the gravitational constant normalized to be 1. Equation pPM2p amounts 
to prescibing the flow to be incompressible, and this is coupled with the 
assumption that there is no flux across the boundary dQ, i.e. v ■ ly = on 

dn. 

As usual, a weak solution to the system pPMip - pPM3p with initial data 
Po € L°°{il.) is defined as a pair {p,v) with 

p e L~(!^ X (0,r)), V e L°°{0,T;L'^{n)), 
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such that for all (j) G C^(M2 x M) we have 




(1) 



(2) 



and for ah V G C^{^) 




(3) 



We remark explicitly that ([2]) includes the no-flux boundary condition for v 
whereas in ([3]) the pressure p has been eliminated (observe also that there 
is no boundary condition on p) . 

Recently, D. Cordoba, D. Faraco and F. Gancedo showed in [I], that on 
the 2-dimensional torus the system OPMip - OPMSp admits nontrivial 
weak solutions {p,v) G L°°(T^ x M) with compact support in time. The 
proof is based on the framework developed in [4] for the incompressible 
Euler equations, although there are several places where the authors had 
to modify the arguments. In technical terms, one of the crucial points in 
the general scheme of convex integration is to show that the relaxation with 
respect to the wave cone of a suitably defined constitutive set, the A-convex 
hull, contains the zero state in its interior. In [T] it was observed, that due 
to a lack of symmetry induced by the direction of gravity, this condition 
seems to fail for IPM; instead, a systematic method for obtaining a suitably 
modified constitutive set was introduced, based on so-called degenerate r4 
configurations. The advantage of the method used in [1] is that it is rather 
robust, and can be used in situations where an explicit calculation of is 
out of reach due to the high complexity (see also \10\ I12j). Indeed, the same 
technique has been recently applied successfully to a large class of active 
scalar equations by R. Shvydkoy [T8] . 

On the other hand, there are certain advantages to obtaining an explicit 
formula for the A-convex hull rather than just showing that a fixed state 
is in the interior. The explicit formula allows one to identify "compatible 
boundary and initial conditions", for which the construction works. For 
the incompressible Euler equations, such initial conditions were called "wild 
initial data" in [5]. For IPM a specific example of independent interest is 
related to the so-called Muskat problem, or 2-phase Hele-Shaw flow. Recall 
that pPMip - OPMSp can be used to model the flow of two immiscible fluids of 
different densities in a porous medium, or, equivalently, in a Hele-Shaw cell 
|17j . If initially the two fluids form a horizontal interface, with the heavier 
fluid on top, it is well known that the initial value problem is ill-posed in 
classical function spaces [22^ [19} [2]. Although some explicit solutions are 
known 0, there is no general existence theory, neither for the evolution 
problem for the interface, nor for weak solutions of IPM. As an application 
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of the explicit formula for the A-convex hull, in this paper we construct weak 
solutions to IPM with initial data given by the unstable interface. 

After a normalization we may assume that the density p{x,t), indicating 
whether the pores at time t near location x G are filled with the lighter 
or the heavier fluid, takes the values ±1. Hence, from now on we will 
concentrate on constructing weak solutions in the sense ([I])- ([3]) which satisfy 

\p{x,t)\ = l for a.e.(x,t) G X (0,r). (IPM4) 

Note that the solutions constructed in [1] also satisfy pPM4p . 

Our weak solutions are in good agreement and show interesting connec- 
tions to predictions concerning the coarse-grained density and the growth 
of the mixing zone made in |141 115j. In |14j F. Otto introduced a relaxation 
approach for pPMip - OFMSp . based on a gradient flow formulation of IPM 
and using ideas from mass transport. It was shown that under certain as- 
sumptions there exists a unique "relaxed" solution p, representing a kind 
of coarse-grained density. Moreover, Otto showed in \15\ Remark 2.1] that, 
in general, the mixing zone (where the coarse-grained density p is strictly 
between ±1) grows linearly in time, with the possible exception of a small 
set of volume fraction 0(t~^/^). 

In Section |4] we show that the relaxed solution p from [T3] is very closely 
related to the concept of subsolution, introduced in [4] for the incompressible 
Euler equations, and in particular we construct weak solutions pk such that 
Pk p. The interpretation is that p is the coarse-grained density obtained 
from yOfc. It is interesting to note in this connection, that, although weak 
solutions are clearly not unique, there is a way to identify a selection criterion 
among subsolutions which leads to uniqueness. 

Moreover, we show that if the coarse-grained density is independent of 
the horizontal direction, the linear growth estimate of [15] is sharp, in the 
sense that there is no exceptional set. As a consequence, we can interpret 
the uniqueness result of Otto as selecting the subsolution with "maximal 
mixing" . In this light it is of interest to note that the analogous criterion 
for the incompressible Euler equations would be "maximally dissipating" 

The paper is organized as follows. In Section [2] we recall from [1] how to 
reformulate (jlPMip - OPMSh as a differential inclusion, and then we calculate 
explicitly the relaxation, more precisely the A-convex hull of the constitutive 
set. These computations form the main contribution of this paper. If one 
is only interested in weak solutions as defined in this introduction (where 
V can be unbounded), the "simpler" computations in Propositions 12.3] and 
13. 11 suffice. However, for completeness we include the computations that are 
required for bounded velocity v in Propositions 12.51 and 13.31 We also show 
that the A-convex hulls obtained are stable under weak convergence, in the 
spirit of compensated compactness [211 E] • 

Then, in Section [3] we show how the explicit form of the A-convex hull 
can be used in conjunction with the Baire category method to obtain weak 
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solutions. For the convenience of the reader we include the details of the 
Baire category method in the Appendix. Finally, in Section |4] we use the 
framework to construct weak solutions to the unstable interface problem. 



2. The relaxation of IPM 

We start by setting 

u:=2v + {0,p). (4) 
Then pPMip - pPMSp can be rewritten as 

dtp + divm = 0, 
div(^x- (0,p)) =0, (5) 
curl(t/ + (0,p)) = 0, 

coupled with 

m = \{pu - (0, 1)), 
|P| = 1 



(6) 



for almost every {x,t). As in [T], we interpret dS])-® as a differential inclu- 
sion: the state variable 



{p,u,m) : X (0,r) ^ M X 



1,2 



is subject to a linear system of conservation laws ([5]), and should take values 
in a constitutive set determined by ([6]). 

The reason for introducing u is purely cosmetic, it makes ([5]) symmetric. 
It also helps in simplifying some of the calculations below. Observe that ([5]) 
can be easily written as a genuine differential inclusion. Namely, we see that 

p - U2 Ui 
Ui P + U2 

is curl -free, hence (locally) equal to for some function cf). Then p = 
and diY {^dtS/ (j) + m) = 0, hence m = —^dt^cj) + V'^ip for some V- In 
particular we deduce that for any {x,t) -periodic solution 

{p, n, m) : ^ M X ]R2 X m2 
of dS]) , there exist periodic functions 

such that 

p=^A(P, 

u=^{2di2(l),d22<p-dn4>), (7) 
m = -^dtV(j) + V^Tp. 
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The periodic solutions of ([5]) are characterized by the wave cone A. This 
is easy to calculate, since it is clearly independent of m and the other two 
equations have a div-curl structure. Hence 

A = {{p,u,m) : IpP = \u\^}. 

We conclude the existence of localized plane-waves, an observation which 
has already been made in [1]: 

Lemma 2.1. Let z = {p,u,m) G A. There exists a sequence 

Wj = {pj,Uj,mj) G {Bi{0) X (-1, x x M^) 

solving ([5]) such that 

(1) dist {wj, [—z, z]) — 7> uniformly, 

(2) Wj ^0 weakly in L? , 

(3) I I\wj\'^dxdt > C\z\^. 

Next, we calculate the A-convex hull of the set in state-space defined by 
([6]). Recall (e.g. from [10], see also [Ml E]) that for a closed set K the 
A-convex hull is defined to be the largest closed set which cannot be 
separated from K: a state z does not belong to if there exists a function 
/ which is A-convex in the sense that 

s I— 7> f{wo + sw) is convex for all u; G A 

such that f < on K and f{z) > 0. It follows immediately from this 
definition that 

zi,Z2 G K with — Z2 € A ^> [^i, 2^2] C K^. 

For further properties of A-convex hulls and functions, see [9]. 
Now, let 

K = {{p,u,m) : \p\ = 1, m = ^pu}. 

Observe that ([6]) is equivalent to 

{p,u,m+ 1(0,1)) e K a.e.. (8) 

We note also that in the absence of spatial boundaries, any solution of ([5]) 
with {p, u,m) G K for a.e.(x, t) is also a solution of ([S])-®. In particular, for 
solutions on the torus we can ignore the additional constant vector ^(0, 1). 

Lemma 2.2. The function 

f{p,u,m) := \m-^pu\ + ^{p'^ + \uf) (9) 

is convex. 

Proof. A short calculation and the triangle inequality shows that, for any 
(p, u, m) and (p, n, m) 

f{p + pt,u + tu, m + tfh) > f{p, u, m) + ct + —{p^ + |up — 2\pu\), 
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where 

c = -{p + u-u) - \fh- \[pu + pu)\. 

The convexity of 

^ ^ fiP + pt,u + tu, m + tfh), 
and hence of / : M x x r2 ^ M follows. □ 



Proposition 2.3. 

=|(/9,'U,m) : \p\ <l,\m- lpu\ < ^{l-p"^)^. 

Proof. Let \pq\ < 1 and uq gM?. Fix e £M? with |e| = 1 and define 

til = uo + (1 - po)e pi = l (10) 

U2 = uo - {I + po)e p2 = -l. (11) 

Then it is easy to see that for any mi , 777-2 

(pi - P2,ui - U2,mi - 7772) G A, 

hence in particular by setting 7r7j = ^piUi we obtain that, with 

7770 : = 1(1 + Po)"ii + ^(1 - Po)m2 
= ^ouo + ^(1 - Po)e, 
the state {pQ,uo,mo) is contained in K^. Thus any {p,u,m) with 

\p\ < 1, 1 777 — i/97t| = i(l — /7^) 

is contained in K^. Moreover, since for any ?f7 € we have (0,0,777.) G A, 
we deduce that in fact any (p, u, m) with 

\p\ < 1, 1 777 — \pu\ < i(l — p^) 

is contained in K^. 

To see that this set is in fact the whole hull, we note that 

Kci{z: g{z) < 1/2}, 

where 

g{p, u, 777) := f{p, u, 777) + -(p^ - \u\^) (12) 

and / is the convex function in ([9]). On the other hand observe that p^ — 
I Tip is A-convex (in fact A-affine). Therefore g is A-convex, and therefore 
necessarily C {z : f{z) < 1/2}. This completes the proof. □ 

Remark 1. Recalling (jS]), we see that if fC denotes the constitutive set given 
by dSD, then (0,0,0) € d)C^. This was the observation made in [U Remark 
4.1]. 
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Proposition 2.4. The set characterizes precisely the relaxation of the 
IPM system. In other words it is stable under weak convergence in the 
following sense: let {pk,Uk,mk) G x L'^ x be a sequence of weak 
solutions to the linear system ([5]) such that 

{Pk,Uk,mk) G for a.e. {x,t), 

and assume {pk,Uk,mi^) {p,u,m) in L^. Then the limit also satisfies 

{p,u,m) S for a.e. {x,t). 

Moreover, is the smallest set containing K with this property. 

Proof. The crucial information is that — |ttp is continuous with respect 
to weak convergence, as a consequence of the div-curl lemma [I3l EH [6] . 
Therefore, if g{pk,Uk,mk) < a.e. (where g is the function from (jl2p ). then 
also g{p,u,m) < a.e. □ 

Although working with the set K would suffice to construct weak solutions 
to ()IPM1|) - ()IPM4|) . there is a slight drawback in that K is not bounded. 
More specifically, the solutions constructed will have u € but not in 
L°°. For the details see the next section. To remedy this problem we next 
consider a compact subset of K. Here the calculation of the A-convex hull 
requires a bit more work. 

Fix 7 > 1 and let 

Kj = {{p,u,m) : \p\ = 1, m = \pu, \u\ < 7}. (13) 

Proposition 2.5. For the set above with 7 > 1, (Kj)^ is given by the 
set of inequalities 

\P\ < 1, (14) 

|n|2<72-(l-p2), (15) 

< 1(1-/52), (16) 



lu\<^l-p), (17) 



|m+in| < |(l+/>). (18) 
Proof. Step 1. Observe that (jl4p . ()17p and (jlSp define convex sets whereas 



(jlSp and (jl6p define sublevel-sets of A-convex functions. Since states in K.y 
satisfy all 5 inequalities, it follows that {K-y)^ is certainly contained in the 
set defined by (fTip - lfTsp . 

Conversely, let be the set of all states (p, u, m) with all inequalities 
([IlD-dlHI) strict. We need to show that C (K-y)^. To this end it suffices 
to show that 

for all z G dU^ \ K-y, there exists z G A \ {0} such that 

z±z£ dUy. ^^^^ 
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Indeed, from this it follows that extr C K^, so that C (^^7)"^ (c.f. the 
Krein-Milman type theorem in the context of A-convexity [H Lemma 4.16]). 
Thus, let z = (p, u, m) G dUy \ K^. 

Step 2. We first note that if (jl6p is an equality, i.e. 



m 



1 1, 2n 

-pu + -{l - p)e 



for some e € 5^, then z lies on a A-segment connecting zi,Z2 € K^. More 
precisely, if 



(l,n+(l-p)e,i(n+(l-p)e)), 
(l + />)e,i(n- (l + p)e)) 



Zl 

then z = i(l + p)zi + ^(1 — /3)z2, as in Proposition! 

Step 3. Next, by elementary computations we verify the following iden- 
tity: 



1 + P 



+ 
1 r 

-u r 



[7^ 



+ 



1 



m H — n 
' 2 ' 



(20) 



Observe that the terms in square brackets are all non-negative in U^. In 
particular, if (jl6p is strict, then either (jl7p or (jlSp has to be strict as well. 
On the other hand, if all three (jl6p - (jl8p are strict inequalities, then it is 
easy to verify that, with z = (1, e, 0), 

Z + tZ ^ U-y 

for all t with \t\ sufficiently small, provided we choose e G 5^ so that p = u-e. 

Step 4. The remaining case is when ()16p and one of (jl7p or (jlSp is a 
strict inequality and the other an equality. By symmetry, we may assume 
that (fTHI) is strict, and 



\m 



i(l-p)- 



Choose e G so that p = u • e, set 



and zt = {pt, ut,mt) := z + tz. Then for all t 

kP=7^-(l-p?) 
\mt - ^ut\ = ^(1 - Pt). 



Moreover, by continuity, (jl6p and (jlSp will continue to hold for |t| small 
enough. This implies (jl9p and thus the proof is completed. 

□ 
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Remark 2. The condition 7 > 1 is sharp. Indeed, if ^ < 1, then it is not 
difficult to see (using standard technology on gradient Young measures) that 
approximate solutions to the corresponding inclusion are compact. In fact, 
ifQ, = T'^, by just looking at pPM2p - pPM3p we deduce easily that any weak 
solution {p,u,m) of ([5])-([6]) with \\u\\l°° < 1 is constant (see for instance 

mm)- 

Remark 3. The computations in Proposition \2.3\ and \2.5\ do not depend on 
the vectors u, m in state-space to he 2- dimensional. Therefore, the same 
formulae continue to hold for the relaxation of the IPM equation in n space 
dimensions with any n > 2. 

3. Construction of Weak Solutions 

There are several ways of constructing weak solutions to differential inclu- 
sions, depending on the particular problems at hand. For the system ([S])-® 
the relaxation is sufficiently large so that a relatively simple iteration proce- 
dure, involving single localized plane- waves, suffices. As always, the crucial 
ingredient is to show that there exists an open set U where states are stable 
only near K (c.f. [U Definition 3.16]). This is expressed by conditions (Hl)- 
(H2) in the Appendix. The methodology of how to pass from this property 
to weak solutions via the Baire category theorem is well known [H [H lU [5] , 
but for the convenience of the reader we present the details in Theorem 15. II 
in the Appendix. 

Recall the definition of K and from Section [21 and let 

U = int K^, = int {K^)^, (21) 

so that 

U =i^{p,u,m) : \p\ < 1, |m - ^pu\ < ^(1 - p^)|, (22) 

and is analogously given by the inequalities (fn|) - (fT8]) all being strict. 

Proposition 3.1. There exists a constant c > with the following property. 
For any {p,u,m) € U there exists {p,u,m) € A with -|- |?fip = 1 such 
that 

{p, u, m) + t{p, u,rh) £ U for \t\ < c(l — p^) 

Proof. From the explicit formulas (jlOp - (jlip it follows that for any {p,u,m) 
with \p\ < 1 and \m — ^pu\ = ^(1 — p^) there exists {p, u, fh) G A with p = I 
such that 

{p, u, m) -|- t{p, u, fh) G for \t\ < co(l — p^) 

for some cq > 0. From this we deduce by continuity the claim (with < c < 
Co) in the case where 

i(l-p2) < \m-^pu\ < i(l-p2). 

In the remaining case, where \m — ^pu\ < \{1 — p^), we can take (p, U, fh) = 
(0, 0, fh) with fh parallel to m — ^pu. □ 
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As a first application we obtain the following variant of [H Theorem 5.2]: 

Theorem 3.2. There exist infinitely many periodic weak solutions to pPMip - 
PPM.SII with 

p e L°°(T2 X M), e L°°(M;L2(t2)), 

such that 

'l a.e. {x,t) G T2 X (0,r), 
fort <0 ort>T. 



\pix,t)\ 



Proof. We construct solutions (p, u, m) of ^ such that 

{p,u,m) G K a.e. {x,t). 

For any such solution = 1 and m = ^pu almost everywhere by definition 
of K, and therefore, on the torus satisfies 

dtp+^dw {pu- {0,1)). 

Therefore, recalling that v = ^{u — (0,p)) we deduce that {p,v) is a weak 
solution of (|IPMip - (|IPM.Sp . 

Next, define the space of subsolutions as follows. Let 

^ = X M, 

^ = X (o,r), 

and 

= |z = ip,u,m) G C^C^) : (P holds and z{x,t) G U V(x,t) G '^}. 

Note that (0, 0, 0) G U, hence Xq is nonempty. Any (p, n, ?7t,) G Xq satisfies 
IpI < 1. Therefore, whenever {p,u,m) G Xq, we have 

\\p\\l^lI{&) < 1) 
hence, using standard elliptic estimates and ([5|), 

for some constant C. Together with the definition of U this implies that 
Xq is a bounded subset of L^j(^) and in particular satisfies condition 
(H3) in the Appendix. Moreover, Lemma 12.11 and Proposition 13.11 imply 
that (H1)-(H2) in the Appendix are satisfied, and consequently Theorem 
2] follows from applying Theorem 15.11 to Xq. 

□ 



Proposition 3.3. Let 7 > 2. If K.y and U-y are as defined in (jlSp and ()2ip . 

then for any z G CZ-y with dist {z, K^) > q > there exists z G A with \z\ = 1 
such that 

z + tz G whenever \t\ < f3, 
where /3 > depends only on a and 7. 
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Proof. Since dU-y is compact and locally the graph of a Lipschitz function, 
it suffices to prove the following quantitative version of (fT9|) : 

for all zq G dU-y \ K there exists e > 0, r > so that 

for ah z G ^^(zo) n dU^ there exists z G A with (23) 

\z\ = r and z it z G dU-y. 

Let Zq = {po, uq, ?no) G dUy \ Ky. 

If |?no — ^PquqI < |(1 — Pq), then using the argument from Step 3 in the 
proof of Proposition 12.51 we may assume that |mo + ^uq] < ^(1 + po). But 
then the A-direction in Step 4 works uniformly for a whole neighbourhood 
of {po,uo,mo). 

Conversely, assume that |mo — ^PoUo\ = ^(1 — Po), so that 
mo = Ipouo + 1(1 - pl)eo with |eo| = 1. 

We know from Step 2 of the proof of Proposition 12.51 that there is a A- 
segment in dU^ of length min(l — po,l + po) centered at {po,uo,mo), but 
now we need to show that the length of the segment can be chosen uniformly 
large for a whole neighbourhood. In other words we need to show that this 
A-segment is stable. 

To this end consider z = {p, u, m) G dU-y nearby, so that in particular 

m = \pu + i(l — p^)e + fh with |e| = 1, 

with I/O— Pol, |tt—uo|, |e—eo|, |?Ti| small. We may also assume that |m—|pu| < 
i(l-p2). Let 

z = [l, e, \e - Y^(m - \u)^ 

as in Step 4 of the proof of Proposition 12.51 and consider zt = {pt,ut,mt) = 
m + tfh. By continuity we can choose e G S""*^ in such a way that 

(p, u, \pu + i(l - p^)e) G dU. (24) 

Since 7 > 2, such e G S*^ exists for any p, u. Moreover, since 7 > 2, the 
circles defined by (I17p - (ll8p are strictly larger than the one defined by (I16p . 
It follows that either (fM|l is satisfied by all e G S^, or the circles defined by 
(|17p -(fT8]) intersect the circle defined by (116p transversely. In particular we 
may choose e and consequently rh so that 

— e • m > 5|m|, (25) 

where (5 > depends only on 7 > 2. 

Now, since e satisfies (f2ll) . {pt,ut) will satisfy (fT5]l for all t G (— (1+p), 1 — 
p) (c.f. Step 2 in the proof of Proposition 12. 5p . Thus, using (i20]l we see that 
Zt remains inside dU^ as long as (I16p remains a strict inequality. Now, a 
short calculation shows that 

1 In 2\ 1 ~ - 

mt - 2PtUt = 2^(1 - Pt)+ -1 _ 
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In particular |mf — ^ptut\ = ^(1 — Pt) if either = 1 or pt satisfies 

(1 + pt){l — p)e ■ fh = — |mp. 

Using (f25]l we deduce that whenever pol, |u— uqI, |e— eo|, |?7t,| is sufficiently 
small, Zt G dU^ for all t with \t\ < l/2min(l — pq,1 + po). This concludes 
the proof. 

□ 

Using the same argument as in the proof of Theorem 13.21 but replacing 
K,U by K^,U^ with 7 > 2, we obtain the following strengthening (corre- 
sponding to Theorem 5.2 of [T]): 

Theorem 3.4. There exist infinitely many periodic weak solutions to pPMl^ - 
PPM3P with p,v e L°°, \p{x,t)\ = 1 for a.e. {x,t) G x [0,r] and compact 
support in time. 

4. Evolution of Microstructure 

In this section we show how to construct solutions to (|IPMip - (|IPM4p in 
the spatial domain 

n := (-1,1) X (-1,1) 

which exhibit the type of mixing behaviour that is expected ([21 HM I15j). 
when one starts with a horizontal interface with the heavier fluid on top. 
Thus, let 

Po{x) = < (26) 

and define subsolutions as follows. 
Let 

& = nx (0,T). 

We consider triples {p,v,m) G L°°(^) such that for all G C^ip. x [0,T)) 

/ I dt4>p + V4>-mdxdt+ [ (l){x,0)po{x) dx = 0, (27) 
Jo Jn Jn 

vV'ipdx = y^p£C°°{^, (28) 



/ 

Jn 



I 

Jn 



n 

{v + (0, p)) ■ V^^p dx = G C^{n), (29) 



and 



\m-pv + i(0, 1 - p^)| < ^(1 - p^) in (30) 
We assume that there exists an open subset ^ C such that 

IpI = 1 a.e. ^\ (31) 

{p,v,m) is continuous in (32) 

\m - pv+l{0, 1 - p2)| < i(l - p^) in ^. (33) 
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Definition 4.1. Let us call any {p,v,m) € L°°{&) satisfying (|27|l - (i33]) an 

admissible suhsolution and the corresponding subset ^ C ^ the mixing zone. 

Next, given an admissible subsolution (j),v,rn), let 

Xq = |(p,n,m) e L°°{&) : {p,v,m) = {p,v,fn) a.e. ^ \ 

and satisfies ([32])- ([33]) in 

As in the proof of Theorem 13.21 the space Xq is bounded in L^(i^), hence 
weak convergence is metrizable on Xq. Let X be the closure of Xq with 
respect to the induced metric. 

Theorem 4.2. There exists a residual set in X consisting of weak solutions 
{p,v) to pPMip - ()IPM4p inQx {0,T) with initial condition pQ. 

Proof. Recall from that v = ^{u — (0, p)), so that (j33p is equivalent to 

{p,u,m + {0,^)) eU for all G 

where U is defined in (j2ip . Therefore Xq satisfies condition (H3) in the 
Appendix. Moreover, (Hi) and (H2) are satisfied by Lemma [2. II and Propo- 
sition [XTl Consequently Theorem 14.21 follows from Theorem 15.11 □ 

Remark 4. If we replace (130p and (j33p by the full set of inequalities (jl4p - 
(jlSp in Proposition \2.5\ for some 7 > 2, and replace U by in the proof 
above, then Theorem remains true with the added information that the 
weak solutions satisfy v € L°°{Sl). 

Remark 5. As a consequence of the residuality in X we obtain, given an 
admissible subsolution (j),v,rn), the existence of a sequence of weak solutions 
{Pk,Vk) of (|IPM1P - (I1PM4P such that 

Pk='pa.e. inSJ\^ and pk -^'p in L°°{^) 

as k 00. In other words p represents a kind of coarse-grained density. 
This justifies calling ^ the mixing zone. 



To conclude this section we exhibit examples of (nontrivial) admissible 
subsolutions. In particular we set 

u = 0, m = (0, m2), 

and assume that (p, m) are just functions of t and the "height" X2, i.e. 

p = p{x2,t), m2 = m2{x2,t). 
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Then (f28]l and (|29]) are automatically satisfied, whereas (f27|) can be written 

as 

dt-p + 92^2 = in (-1, 1) X (0, T), (34) 
rn2 = for X2 = ±1, (35) 

[-1 X2<0 , 
p = < for t = 0, 36) 

^ [ + 1 X2>0 ^ ' 

interpreted in the weak formulation. Furthermore, (j30p becomes 

- (1 -P^) < "i2 < in (-1,1) X (0,r). (37) 

Note also in connection with Remark S] that in this case ()15l ). ()17p and (llSp 
follow automatically from (|37p provided 7 > 3. 

An obvious way to construct admissible subsolutions is then to prescribe 

m2 = -a(l - p^) 

for a G (0, 1). The resulting equation from (I34p is essentially the inviscid 
Burgers' equation, with the (unique) entropy solution given by 

{—1 X2 < —2at, 
1^2 1 < 2at, (38) 
+1 X2 > 2at 

for times t < l/(2a). The corresponding mixing zone is 

^ = {{x,t) G ^ : \x2\ < at}. 

Since at t = l/(2a) we have |p| < 1 for all X2 € (—1, 1), the functions p, m 
can easily be extended continuously to later times so that (|34p - (|37p continue 
to hold. 

It is interesting to note that in the borderline case a = 1 the subsolution 
p is precisely the unique solution of the relaxation approach to the problem 
(|IPMip - pPM4p introduced by F. Otto in [H]. Thus, although weak solu- 
tions are clearly not unique, there seems to be a way to recover uniqueness at 
least for subsolutions. We plan to explore further this connection elsewhere. 
Here we contend ourselves with showing that among all subsolutions, for 
which p is a function of t and the vertical direction X2 only, the case a = 1 
in (j38p corresponds to the one with "maximal mixing". This gives a new 
interpretation to the results in 



Proposition 4.3. Let {'p,v,rn) be an admissible subsolution to the unstable 
initial condition (j26p such that p = p(x2,i)- Then the mixing zone ^ is 
contained in {{x,t) : \x2\ < 2t}. 

Proof. We start by observing that since dip = 0, (|28p - (|29p imply that v = 0. 
Next, let 

^0 s < 0, 



s s > 0, 
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and consider the test function (j){x, t) = x{x2 — 2t) in (f27|l (it is easy to see 
by approximation that this is a valid test function). We obtain 



m2 — 2pdxdt+ / X2Po{x)dx = 0. 



'{{x,t)€&':X2>2t} J{xGQ:X2>0} 

Therefore 

(m2 -2p + 2) dxdt = 0. 



'{{x,t)£S>:X2>2t} 

On the other hand ()30p imphes 

^2 >-(!-/). (39) 

We deduce that 

p = 1 a.e. {x,t) e {x2 > 2t}. 
Analogously we also find 

p = — 1 a.e. {x,t) G {x2 < —2t}. 



□ 



5. Appendix 

We consider general systems in a domain & C M"^ of the form 



AAz = in ^ (40) 



i=l 

z{y) G K a.e. y G & (41) 

where 

is the unknown state variable, Ai are constant mx N matrices, and K C 
is a closed set. We make the following assumptions. 



(HI) The Wave Cone: There exists a closed cone A C and a constant 
C such that for all z G A there exists a sequence wj G C^(i?i(0); M^) 
solving (|40p such that 

(1) dist {wj, [—z, z]) — > uniformly, 

(2) Wj weakly in L^, 

(3) J\wj\'^dy > C\zf. 

(H2) The A-convex hull: There exists an open set \J C with \J n 
K = and such that for all z G U with dist (z, K) > a > there exists 
z G A n S^^^ such that 

z + tzeU ioT: ah |t| < /3, (42) 

where /3 = /3(a) > 0. 
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(H3) Subsolutions: Xq is a nonempty bounded subset of L^(^) consisting 
of functions which are " perturbable" in an open subdomain ^ d . This 
means that any z G Xq is continuous on ^ with 

z(yy) eU for y e ^, (43) 

and moreover, if z G Xq and w G Cc{^) such that w solves (jlO|) and 
{z + w){y) G [/ for aU y G then z + w e Xq. 

Finahy, let X be the closure of Xq with respect to the weak topology. 
Since Xq is bounded, the topology of weak convergence is metrizable on 
X, making it into a complete metric space. Denote its metric by dx{-, •)• 

Theorem 5.1. Assuming (H1)-(H2), the set 

{z e X : z{y) G K a.e. y G '^j 

is residual in X . 

The proof relies on the following lemma, where, for notational convenience 
we set 

F{z) := min(l,dist iz,K)). 

Lemma 5.2. Let z G Xq with F{z{y))dy > e > 0. For all rj > there 
exists z G Xq with dx{z, z) < rj and Ja^^ \z — z\'^dy > 6, where 6 = 6{e) > 0. 

Proof of Lemma \5.2[ Since z G Xq is continuous and z{y) ^ K on for 
any yo G there exists tq = ro(yo) > such that 

\F{z{yQ)) < F{z{y)) < 2 F(z(yo)) (44) 

for all y G Brgiyo) C Then, by a simple domain exhaustion argument, 
we find disjoint balls Bi := Bri{yi) C ^ for i = 1 . . .k such that 

Im, (45) 

i / Fiz{y))dy. (46) 

Next, observe that (H2) implies the existence of a continuous function 
4> : [0, 1] [0, 1] such that (p{0) = 0, (j){t) > for i > and such that for 
any z (zU there exists z G A n S^~^ with 

z + tz whenever |t| < 4i[F{z{y))). (47) 

By considering its convexification if necessary, we may assume without loss 
of generality that cj) is convex and monotone increasing. Using (|44p . ()46p 



i=l 

J2 I nz{y))dy > 

JB,, 
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and the convexity of (j) we obtain 



m 

\Bi] 
m 



Moreover, using (jl7|) . (HI) and a simple rescaling, there exist Zj € C^{Bi; 
such that 

(1) z{y) + Zi{y) G U for all y, 

(2) dizi,0) <ri/k, 

(3) Jj,Jzi\^dy > C(P{Fiz{yi)))\B,\, 

where C is the constant in (HI). Therefore, 



w 



:=^5i eCe(^) 



and z{y) + u;(y) € [/ for all y € hence z := z + w € Xq hy (H3). 
Moreover, z £ Xq satisfies 

dxiz,z) < y^^dxizi,0) < 7] 

i 

i 

This concludes the proof. □ 

Proof of Theorem \5.1[ First of all the functional I{z) = J^^ \z{y)\'^ dy is a 
Baire-1 function on X. Indeed, observe that 



Ij{z) := / \z* pj{y)\ dy, 

J{ye^^:dist iy,d'?/)>l/j} 

where pj G C^{Bi/j{0)) is the usual mollifier sequence, is continuous as a 
map X ^ M, and moreover 

Ij{z) — I{z) as j oo. 

Therefore, by the Baire category theorem the set 

Y := {z X : I is continuous at z} 

is residual in X. We claim that z implies J^^ F{z{y))dy = 0. 

If not, let e := F{z{y))dy > for some z £ Y, and let Zj € Xq be 

a sequence such that Zj z in X. Since I is continuous at z, it follows 
that Zj — >■ z strongly in L^(^), and in particular we may assume that 
J^Fiz^iy))dy>e/2. 
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Then, by applying Lemma 15.21 to each Zj, we obtain a new sequence zj £ 

Xq such that Zj z in X (and hence strongly in L^), but \zj — Zj\^dy > 
5 > 0, where 6 only depends on e. This contradicts the strong convergence 
of both sequences to z. 

□ 
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